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ABSTRACT 

Stars close to the Eddington luminosity can have large low-density inflated envelopes. We show that the rise times of shock breakout 
signals from supernovae can be extended significantly if supernova progenitors have an inflated stellar envelope. If the shock breakout 
occurs in inflated envelopes, the shock breakout signals diffuse in them, and their rise time can be significantly extended. Then, the 
rise times of the shock breakout signals are dominated by the diffusion time in the inflated envelope rather than the light-crossing 
time of the progenitors. We show that our inflated Wolf-Rayet star models whose radii are on the order of the solar radius can have 
shock breakout signals that are longer than ~ 100 sec. The existence of inflated envelopes in Wolf-Rayet supernova progenitors 
may be related to the mysterious long shock breakout signal observed in Type Ib SN 2008D. Extended shock breakout signals may 
provide evidence for the existence of inflated stellar envelopes and can be used to constrain the physical properties of these enigmatic 
structures. 

Key words, stars: evolution - stars: massive - supernovae: general - supernovae: individual (SN 2008D) 


1. Introduction 

The first electromagnetic signals from a supernova (SN) 
are emitted when the shock wave that travels thro ugh the 
stellar interior reaches the stellar surface (e.g., Cplgafi 


19741: Klein & Chevalied 1978t Ensman & Burrows 19921 


type and the early light curve modeling p oint to a WR progeni- 
tor with a radius smaller than 10 fe.g.. lCh evalier & Eryssonl 


Colgatd 

Il992t 


,_, _ „,la _ 

2008L iRabin^ & Waxmanll201 iL see also iDessart et al.l 1201 iL 
Bersten et al.l 2013 ). The observed long shock breakout duration 


Matzner & McK^ 1999t iTominaga et al. T 2 OI it iTolstov et afl 

2 OI 3 I) . Photons trapped in the shock wave are suddenly released 
near the surface and start to be observed at this moment. This 
sudden release of photons in the shock wave near the stellar sur¬ 
face is called shock breakout. 

The shock breakout signals are strongly affected by the prop¬ 
erties of the SN progenitor and can be used to estimate them. 
Especially the duration of the shock breakout signals is assumed 
to approximately correspond to the light-crossing time fie of the 
progenitor, that is, fie - /?*/c, where /?* is the progenitor radius 
and c is the speed of light. If a Wolf-Rayet (WR) star explodes, 
its shock breakout signal is expected to have a duration of about 
1-10 sec beca use typical radii of WR stars are on the order of 
1 - lORo (e.g.. ICrowtheill2nn7l lYoon 1 Oh . 

The shock breakout durations are, however, suggested to 
be extended for several reasons. For example, if the circum- 
stellar medium (CSM) density is high enough to be optically 
thick, shock breakout durations can be exte nded as a result 
of ph o ton diffusion in t he den se CSM (e.g., iBalberg & LoebI 
I 2 OI II: ISvirski & Nakad l2014allbll . Aspherical photospheres 
may also cause an extension of the shock br eakout durations 
dSuzuki & Shigevamall2010l:ICouch et al.ll20li]) . In addition, the 
duration of the shock breakout signals from Type Ib SN 2008D 
was observed to last for ~ 300 sec with a rise time of ~ 60 sec 
dSoderberg et al.l2008l:lModiaz et al.ll2009l) , indicating a progen¬ 
itor radius of R* ficC ~ 130 Rq. The indicated progenitor ra¬ 
dius is too large to correspond to a WR star, but its SN spectral 


may also indicate the existence of mechanisms to extend shock 
breakout signals. 

In this Letter, we argue that shock breakout rise times can be 
extended significantly if SN progenitors have inflated envelopes. 
In that case, the shock breakout signals diffuse in the inflated 
envelope, resulting in the extension of the shock breakout rise 
times. The diffusion times fd of canonical WR SN progenitors 
are shorter than the light-crossing time fic (Table [T]l. However, 
we show that the diffusion time can be much longer than the 
light-crossing time in inflated WR SN progenitors even if their 
radii are on the order of the solar radius. 

We first explain the inflation of stellar envelopes and our in¬ 
flated stellar evolution models in Sect. We show the effect 
of the inflated envelopes on the SN shock breakout signals in 
Sect. |3] Our conclusions and their implications are summarized 
in Sect. |4] 

2. Inflation of stellar envelopes 

We first discuss the inflation of stellar envelopes. Inflation 
refers to the extended low-density envelopes in stel lar models 
that r e ach Eddington lumin o sity in their outer laye rs (llshii et al.l 


I 1999 I: iPetrovic et al.l 120061: iGrafener et al 


nomenon was first discussed bv iKatoT i 1985 


20121 


198 


The phe- 
in the con¬ 
text of supermassive stars, considering electron-scattering as 
the only opacity source. The envelope inflation in the models 
discussed below is, how ever, caused by t he Fe-group opacity 
bump at T ~ 200 kK (llshii et al.l 119991: IPetrovic et al.l 1200^ 
IGrafener et al.l 1201^ . When the stellar envelope approaches 
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the Eddington luminosity, it will become convectively unstable 
(lLangerlll997h . The convective luminosity would need to carry 
any super-Eddington flux so that the radiative luminosity does 
not exceed the Eddington luminosity. However, if the envelope 
density is too low, convection cannot become efficient enough 
to transport sufficient energy inside the envelope. Then, the star 
expands its envelope in order to find a structure with a lower 
opacity so that the Eddington luminosity is not exceeded. 

To discuss the actual structure of inflated envelopes in WR 
SN progenitors, we performed stellar evolution calculations 
by using an advanced one-dime nsional hydrodynamic binary 
stelly evolution code (EEC; see iHeger et al]l2000l: lYoon et ^ 
120061 1201 nl: [Rrott et al Jl2nTl]f to model WR stars as nonrotat¬ 
ing hydrogen-free helium stars. The code solves the Lagrangian 
version of the equations for mass, energy, and momentum con¬ 
servation together with convective and radiative energy trans¬ 
port and is explicitly coupled t o a diffusive mixing s cheme and a 


nucleosvnthesis network (seelKozvrevaet al 

20141 for details). 

BEC incorporates the latest input physics ( 

Heser et all 

2000 t 


lYoon et al.l20Ml2010HBrott et alJ201 IbfKohler et aAoiSh and 


is well-suited to investigate st ars evolving close t o the Eddington 
limit with inflated envelopes (iKohler et al.ll2015l ; D. Sanyal et al. 
in preparation). 

In particular, the opacities used to compute the models wer e 
interpolated from the OPAL tables (llglesias & Rogerj Il996ll . 
which contain an opacity enhancement around T ~ 200 kK due 
to Ee-group elements. Convection was tre ated within the frame- 
work of the standard mixing-length theory (iBohm- Vitensel 1 95^ 
with the mixing length set to 1.5 times the lo cal pressure scale 
height . The WR mass-loss recipe proposed bv iNugis & Lamei^ 
(l2000l) was used for our computations. The stellar surface was 
set at the photosphere where the optical depth is 2/3. We inves¬ 
tigate the evolution of helium stars with initial masses of 10 and 
12 Mq in this Letter. We explored two sets of models, one with 
solar metallicity and the other with half-solar metallicity. The 
half-solar metallicity models were computed until the oxygen¬ 
burning stage, while solar metallicity models were only com¬ 
puted until carbon ignition because we encountered numerical 
problems to iteratively converge models with large inflated en - 
velopes that are larger at higher metallicity (llshii et al.lll999l) . 
The inflation in these models at the time of the SN explosion 
could thus be somewhat stronger than we assume here. The evo¬ 
lutionary tracks of the models in the Hertzsprung-Russell (HR) 
diagram are presented in Eig. [T] The final stellar properties are 
summarized in Table [1] 

The helium zero-age main-sequence models of both 10 and 
12 Mq stars are marginally inflated. As a consequence of the ap¬ 
plied mass loss, they lose mass and decrease in luminosity dur¬ 
ing the core helium-burning phase. After helium is exhausted 
in the core, the stars begin to contract and become hotter and 
brighter. The models are hardly inflated during this phase be¬ 
cause the Ee-group opacity bump is only partially contained in¬ 
side the stars. When the helium-shell ignites, the evolutionary 
tracks eventually move t oward cooler temperatur es because of 
the mirror principle (e.g.- lKinnenhahn et al.ll201^ . As the mod¬ 
els become cooler, they become significantly inflated and ex¬ 
hibit a pronounced core-halo structure, as shown in Eig. |2 The 
electron-scattering Eddington factors of our models are shown 
in Table [T] When the full opacity is counted into the Edding¬ 
ton factor, all models reach the critical value of one. The en¬ 
velope inflation is stronger for the solar metallicity models be¬ 
cause the prominence of the Ee-group o pacity bump incre ases 
with metallicity, as fo und in the studies bv llshii et al.l(ll999h and 
iPetrovic et al.l (l2006t) . Eor comparison, we also show a poly¬ 



Fig. 1. Evolutionary tracks of the helium star models in the HR dia¬ 
gram. The initial metallicity, initial mass, and final mass of the models 
are shown. The locations of helium-shell ignition (circles), core car¬ 
bon ignition (squares), and core oxygen ignition (triangles) are marked. 
Lines of the constant radii of 0.5, 1, and 2 Rq are also shown. 



radius (Rq) 


Fig. 2. Final density structure of our stellar models. We also show the 
structure of a star without the envelope inflation, which is a polytiropic 
star of 7 Mq and 2.4 Rq with the polytropic index of 3. The inflated 
stellar models have the extended low-density regions on top of the core 
structure. 


tropic star with M* = 7 Mq and R* =2.4 Rq obtained with the 
polytropic index of 3 which does not have an inflated envelope. 
Figure|2]shows that our stellar models have extended low-density 
layers on top of the core structure. This type of the envelope in¬ 
flation is also found i n the previously re ported models with sim¬ 
ilar final core masses dYoon et al.ll2012h . The evolution and final 
position of our stars in the HR diagram are consistent with their 
models. 


3. Shock breakout in inflated stellar envelopes 

The shock breakout in a SN occurs when photons start to leak 
from the shock wave, that is, when the dynamical timescale fdy 
of the shock propagation in the remaining unshocked stellar en¬ 
velope with radius AR (fjy = AR/vsh) becomes similar to the 
diffusion timescale fd = tAR/c, where Vsh is the shock velocity 
and r is the optical depth in the remaining unshocked stellar en¬ 
velope (e.g.. lWeaveilll97^ . The shock breakout condition can 
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Table 1. Stellar properties and corresponding timescales: Initial metallicity, initial stellar mass, final stellar mass, final surface electron-scattering 
Eddington factor, final mass-loss rate, final ratio of the mass-loss rate to the critical mass-loss rate, final stellar radius, light-crossing time fic = 
distance from the radius at t = 30 to R^,, diffusion time for t = 30 (v^h = 10000 km s“'), distance from the radius at t = 20 to diffusion time 

for T = 20 (Vsh = 15000 km s“'), distance from the radius at t = 10 to f?*, and diffusion time for t = 10 (vsh = 30000 km s“*). 


z 

(Zo) 

Mi 

(Mo) 

M* 

(Mo) 

F. 

M 

(lO-'-Moyr-') 

^ / -^max 

R* 

(Re) 

^Ic 

(sec) 

AR^" 

(Rs) 

(sec) 

AR^“ 

(Re) 

(sec) 

AR'" 

(Re) 

4*' 

(sec) 

1 

10 

6.8 

0.412 

8.1 

0.45 

2.43 

5.67 

1.73 

121 

1.01 

47.0 

0.130 

3.01 

1 

12 

7.8 

0.439 

7.6 

0.57 

2.31 

5.38 

1.71 

119 

1.70 

78.6 

0.372 

8.62 

0.5 

10 

7.4 

0.398 

6.4 

0.066 

1.06 

2.47 

0.189 

13.2 

0.188 

8.72 

0.163 

3.78 

0.5 

12 

8.6 

0.428 

8.3 

0.11 

1.33 

3.10 

0.356 

24.9 

0.351 

16.27 

0.306 

7.09 

- 

- 


- 

- 

- 

2.4 

5.60 

0.0187 

1.31 

0.0170 

0.788 

0.0140 

0.325 


Notes. Polytropic star with the polytropic index of 3. 



Fig. 3. Optical depth t in our models. The optical depth in the 
polytropic star is obtained with the constant opacity of 0.2 cm^ g“'. 


be simply expressed as t c/vsh- For a shock breakout at the 
stellar surface without the inflated envelope, AR is so small that 
fd is smaller than the light-crossing time t\c - R-k/c, making the 
shock breakout duration - /?*/c. For example, in the case of 
the 2.4 Rq polytropic star we show in Fig. |2l the light-crossing 
time is fie - 5.60 sec and the diffusion time is fa - 1.31 sec, 
assuming Vsh = 10000 km s“' (Table[T]i. 

The inflated stellar models are characterized by extended 
low-density envelopes (Fig. |2]l. As a result, AR can be much 
larger than the corresponding values for the uninflated stellar 
models. In Fig. |3l we show the optical depth in our stellar mod¬ 
els. In Table [1] we summarize the light-crossing time fic and the 
diffusion time fa for several shock velocities. In the uninflated 
polytropic model, the shock breakout occurs at the edge of the 
stellar surface (AR — 0.0187 Rq for r = 30), making the diffu¬ 
sion time shorter than the light-crossing time. If the progenitor 
has an inflated envelope, however, the shock breakout can occur 
at the bottom of the inflated envelope, keeping AR large when 
the shock breakout occurs. For example, our solar-metallicity 
10 Mq model has AR - 1.73 Rq for r = 30 because of the enve¬ 
lope inflation. The large AR makes the diffusion time (fa - 121 
sec) much longer than the light-crossing time (fic - 5.67 sec). 
Thus, the timescale of the shock breakout signals within the in¬ 
flated envelopes can be dominated by the diffusion time, and 
their rise time is determined by the diffusion time. The subse¬ 
quent light curve is expected to decline exponentially with an 
e-folding timescale of the diffusion time due to the photon dif¬ 
fusion in the shocked envelope. Overall, the total shock break¬ 


out duration becomes significantly longer than the light-crossing 
time. The dynamical timescales of our inflated envelopes at the 
time of their explosions are ~ 2000 sec, and the inflated en¬ 
velopes can be sustained until the shock passes through. 

We have computed the diffusion times for several shock ve¬ 
locities between 10000 and 30000 km s '. If the shock velocity 
at the shock breakout is too high, the diffusion time becomes 
shorter than the light-crossing time. In our models, the diffusion 
time starts to be shorter than the light-crossing time when the 
shock velocity is higher than ~ 30000 km s * (Table [1]). The 
shock velocity is suggested to reac h above 30000 km s ' at the 
shock breakout in WR stars (e.g.. iNakar & Sarill2010h . How¬ 
ever, previous studies assume that WR stars only have a steeply 
declining envelope as in our polytropic model. A flat inflated 
envelope exists on top of the steeply declining core in the den¬ 
sity structure of our models. The flat density structure deceler¬ 
ates the shock significantly, and the shock velocity at the shock 
breakout in inflated envelopes can be lower than 30000 km s '. 
Because the shock propagation in inflated stellar structure has 
not yet been investigated, hydrodynamical modeling in inflated 
envelopes is required to estimate the shock velocity at breakout 
in this case. 

The sizes of the inflated envelopes in the half-solar- 
metallicity models are smaller than those in the solar-metallicity 
models because the Fe-group opacity bump is weaker with de¬ 
creasing metallicity. However, the diffusion time can still be 
longer than the light-crossing time, and the shock breakout rise 
times can be longer than 10 sec in these WR stars. 

Although higher metallicity models are preferred to have 
large AR and thus long fd in terms of opacity, the inflated 
envelope can disapp ear because of higher mass-loss rates. 
iPetrovic et alJ ( l2006l) showed that the inflated envelope disap¬ 
pears when the mass-loss rate of the star is higher than the criti¬ 
cal mass-loss rate defined as 


o / KjiVl * 

Afmax - ^JrRlpm J — -, ( 1 ) 

V 

where is the lowest density in the infl ated envelope and 
is the radius at the lowest density (see also lGrafener et al]l2012l) . 
We show the ratio of the mass-loss rate in our models to the 
critical mass-loss rate in Table [1] Our solar-metallicity models 
are close to the critical mass-loss rates. Higher-metallicity stars 
are likely to exceed the critical mass-loss rate. 

Although we have focused on WR SN progenitors in this 
Letter, stars do not need to be WR stars to have inflated en¬ 
velopes. When stars are near the Eddington luminosity, they 
can have inflated envelopes even if they are hydrogen-rich. For 
example, luminous blue variables are hydrogen-rich stars that 
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are close to the Eddington luminosity, an d they have recently 
been suggested to be SN progenitors (e.g.. lGal-Yam & LeonardI 
|2009|) . They can also have inflated envelopes when they explode. 

To confirm that some SN progenitors are actually inflated 
stars, observational consequences of the inflated envelopes in 
SN observational properties other than the shock breakout rise 
times are required to be investigated. Because the inflated en¬ 
velopes in our models only have ~ 10 “^ Mq, they are not likely 
to strongly affect later SN observables such as its light curve. 
However, the shock breakout sp ectra or very early SN spectra 
such as that recently obtained by iGal-Yam et al.l (1201 4h may be 
affected by an inflated envelope. It is especially important to 
distinguish the shock breakout extension caused by an inflated 
envelope from that produced by a dense stellar wind because the 
shock breakout signals can be extended by the photon diffusion 
in both cases. Estimates of the progenitor mass-loss rate, for 
instance from multiwavelangth SN observations, may constrain 
the wind density near the stellar photosphere and thereby con¬ 
firm or rule out a shock break-out e xtension caused by winds (see 
Sect. B.ll for the case of SN 2008D) . [iVIaedal (|20 1 3h also discusses 
possible observational consequences of the inflated envelopes. 


3.1. SN2008D 


The extension of shock breakout signals by an inflated envelope 
may explain the long duration of the shock breakout signal ob¬ 
served in SN 2008D. Even if the progenitor of SN 2008D is a 
WR star with R* ~ Rq, our result show that the rise time of 
the shock breakout signals can be ~ 60 sec because of the long 
diffusion time in the inflated envelope below which the shock 
breakout occurs. The shock-breakout luminosity decline after 
the luminosity peak is found to be sim ilar to an exponential de¬ 
cay for a while (ISoderberg et al.ll2008l) . which is consistent with 
the diffusion interpretation. 

The observed high temperatu re of the shock breakout sign als 
in SN 2008D (-0.1-1 keV, e.g.. lModiaz et al.l2009llLil2(X^ is 
sugg ested to indicate t he shock velocity abov e ~ 20000 km s"* 
(e.g.. iKatz et akll^OlOl : [Balberg & Loeb(l201 Ih . The longest dif¬ 
fusion time for Vsh = 20000 km s ' among our models is 19 sec 
from the solar 12 Mq model. Although it is shorter than ob¬ 
served, ~ 60 sec, the diffusion time is still much longer than the 
light-crossing time. The o bserved e-folding time scale in the de¬ 
clining phase (~ 130 sec. ISoderberg et al.ll2008h is also longer 
than the diffusion time in our models, but our models do pre¬ 
dict the overall extension of the shock breakout signals. The 
shock temperature becomes ~ 0.1 keV using v^h = 20000 km s“* 
and our inflated envelope density (iKatz et al.l2010l) . which corre- 


' t al.l2010h. ’ 

sponds to the observational constraints (e.g.. lModiaz et al.l2009l) 
and can explain the observed thermal X-ray luminosities (~ 10"'^ 
erg s“') with a radius of - Rq. More investigations are required 
to see whether the X-ray spectra can be explained by our model. 

The later shock velocity is observed to be 0.25c 
dSoderberg et alj|2008l) . However, the shock can have acceler¬ 
ated after the shock breakout because of the steep density de- 
cline between the s tellar surface and the low density wind (cf., 
iTolstov et alj|2013l) . Thus, the observed velocity may not corre¬ 
spond to the actual shock velocity at the shock breakout. 

The estimated SN ejecta mass of SN 2008D is 3 - 7 Mn 
( Soderberg et alJ 200^ iMazzali et alj|2008t iTanaka et all 1200^ 


iBersten et al. 12013 ). Our pre-SN stellar models (Table [T]i have 

6.3 - 7.8 Mq. Subtracting a remnant mass of about 1.4 Mq, 
the expected ejecta masses from our inflated stellar models are 
consistent with those of SN 2008D. The mass-loss rates of our 


models are also consistent with that estimated fro m radio obser¬ 
vations (7 X 10 ® Mq vr~'. ISoderberg et all2008h . 

The metallicity at the location where SN 2008D appeared is 
estimated to be 0.5 Zq based on th e metallicity gradient in the 
host galaxy dSoderberg et al.ll2^8ll . Our 0.5 Zq stars do not in¬ 
flate enough to explain the long shock breakout rise time of - 60 
sec. Howe ver, Galactic metallic ity gradients have large disper¬ 
sions (e.g.. lPedicelli et al1l2009l) . Since our Zq stars have much 
longer diffusion times, it is possible that the progenitor metallic¬ 
ity of SN 2008D is around Zq or higher, and the progenitor had 
a sufficiently inflated envelope to show the long shock breakout 
rise time with the estimated shock velocity. Recent metallicity 
measu rements at the location of SN 2008D support our argu¬ 
ments dModiaz et al.ll20lTllTh6ne et alJl2009l) . 


4. Conclusions 

We have shown that the rise times of supernova shock break¬ 
out signals can be extended because of inflated stellar envelopes. 
The shock breakout can occur within the low-density inflated 
envelopes in which the shock breakout signals are diffused and 
extended. The long diffusion time in the inflated envelopes 
makes the shock breakout rise times long. The shock breakout 
timescale is then dominated by the diffusion time instead of by 
the light-crossing time. Even if a SN progenitor has a radius on 
the order of the solar radius whose light-crossing time is a few 
seconds, the rise time of the shock breakout signals can be more 
than 100 seconds because of the inflated envelope. The exten¬ 
sion of the shock breakout signals by the inflated envelope may 
explain the mysterious long shock breakout rise time observed 
in Type Ib SN 2008D, although more investigations are required 
to confirm this. 

The actual existence of inflated stellar envelopes is still 
widely debated. While we And the envelope inflation in our stel¬ 
lar evolution models, the envelope inflatio n is not as prominent 
in other models (e.g.. lEkstrom et aO2012l) . If infla ted envelopes 
exist, they are also pulsationally unstable (e.g., iGlatzel et al.l 
1 19931) with unknown consequences. If we can confirm that SN 
shock breakout signals from WR SN progenitors have signatures 
of inflated envelopes, for example, extended shock breakout sig¬ 
nals from SN progenitors with low mass-loss rates, it may indi¬ 
cate the common existence of inflated stellar envelopes in nature. 
The shock breakout signals can be a viable probe of unsolved 
problems in stellar structure as well as a touchstone of the stellar 
evolution theory, which predicts the final status of the massive 
stars. 
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